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I. INTRODUCTION 


This report describes a computer program which analyzes im- 
pedance and pattern characteristics of wire type antennas mounted 
on a satellite. The satellite body is modeled by an enclosing struc- 
ture of flat plates, some of which act as ground planes for the 
wire antennas. The basic principles and techniques for the present 
analysis are similar to those developed previously for the case 
of wire antennas on a flat plate! This work extends the analysis 
of reference 1 to account for a more complex multiplate structure 
and deals with the general case where several antennas can share 
the same ground plane or they can be located on different faces 
of the satellite structure. As a result, to obtain accurate values 
for the antenna impedances, it is essential that all coupling terms 
due to reflected and diffracted fields from nearby panels and wedges 
as well as those due to the presence of other wire antennas must 
be added to the appropriate terms of the open-circuit impedance 
matrix. Efforts have been directed to make the present program 
compatible in programming logic and terminology with existing thin 


wire antenna programs developed by Professor J. Richmond2?3, 


II. INPUT DATA AND SUBROUTINE CWIRE 


Subroutine CWIRE is used to setup the input data for the main 
program. The following parameters must be specified in subroutine 
CWIRE 


AL wire radius in wavelengths 

CMM wire conductivity in megamhos/m 

FMC frequency in MHz 

NPGP number of wire end-points on the ground planes 

NRP number of real points including those on the ground planes 
NSGP number of wire segments touching the ground planes 

NRS number of real wire segments. 


To illustrate the working of the subroutines and main program, 
a satellite structure consists of 10 plates as sketched in Figure 
1 is used in the examples. In the first example, a quarter wave 
monopole is mounted at (0.,-A,0.) on plate 1 of the structure. Sub- 
routine CWIRE for this case is listed in Figure 2. ‘n the second 
example, another quarter wave monopole is located at the center 
of plate 2 and perpendicular to this plate, the corresponding CWIRE 
is listed in Figure 3. In both cases, CWIRE first defines the wire 
geometry by specifying the coordinates XC, YC, ZC of the endpoints 
in meters. For more complex wire structures such as that of helical 
antennas, a subroutine which generates XC, YC, ZC similar to sub- 
routine HELIX of Appendix 2 can be called for this purpose. Do 
loop 702 generates the segment numbers and the list of end-points 
IA, IB for each real wire segment (Figure 2). The array element 
IDPT (1) specifies the plate (number) that acts as the ground plane 
for endpoint I. It then reads in and stores input data describing 
the model of the satellite structure, which include the position 
vectors CR of the corners and the list of corner index NPLC. Figure 


1 shows the carrier index written next to the corner, numbered from 
1 to 12. The plate number is written on each plate and encircled. 
Plate 10 is the bottom plate defined by corners 9, 10, 11 and 12. 
This input data for satellite geometry is labeled file IN] in user 
name 3468N. Type - COPY IN1,3468N TO .IN, user name will make this 
file available for input to subroutine CWIRE. The data file INI 

is listed below. 


The meaning of the above figures is self-explained in the fol- 
lowing tables. 


Table 1 
Corner index x y 
ais) ain (in wavelengths) 
CR(I,1) CR(I,2) CR(I,3) 

1 1.5 -2.0 0.0 
2 1.5 15 0.0 
3 71.0 1.5 0.0 
4 -1.0 -2.0 0.0 
5 208 -3.0 -1.0 
6 2.5 3.0 -1.0 
7 -2.5 3.0 -1.0 
8 -2.5 -3.0 ~1.0 
9 2.5 -3.0 -2.5 
10 2.5 3.0 2.5 
11 -2.5 3.0 2.5 
12 -2.5 -3.0 =2.5 


Each plate is defined by 4 corners Cl, C2, C3, C4. 


Table 2 
Plate Number Corner Index 
Cl C2 C3 C4 
NPLC(J,1) NPLC(J,2) NPLC(J,3) NPLC(J,4) 

1 1 2 3 4 
2 5 6 2 1 
3 2 6 7 3 
4 4 3 7 8 
5 5 1 4 8 
6 9 5 8 12 
7 9 10 6 5 
8 6 10 11 7 
9 8 7 dil 12 
10 9 12 11 10 


The position vectors CR of the corners are converted to metric 
units in Do 18. Subroutine CWIRE then generates the list of unit 
vectors VNP, where VNP (I,3) denotes the unit vector normal to plate 
I. The normal unit vector is defined as a cross-product of two 
edge vectors and points outward, away from the structure. This 
is done in loop Do 12. Note that the corner index NPLC (I,J), J= 
1,2,3,4 has been arranged such that an observer proceeds from one 
corner to the next on the same plate in the direction of increasing 
J will see the normal vector to his left, pointing outward. 


| 


Figure 1. Model of 10-plate structure used in examples. 
Dimensions in wavelengths. 
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Figure 2. Subroutine CWIRE. 
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Figure 2. (Cont'd). 
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Figure 3. Subroutine CWIRE. 
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Figure 3. (Cont'd). 
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IIT. THE MAIN COMPUTER PROGRAM 


The main computer program is listed in Figure 4. This program 
calls subroutine CWIRE for input data, namely the geometry of the 
wire antennas and the multiplate structure. The program next calls 
subroutine CSORT to generate appropriate image points and image 
segments. CSORT also sets up the dipole current modes for the wires 
and their images. Each dipole mode I has segments JA(I) and JB(I), 
terminals at point I,(1), and endpoints 1,(1), 1,(1). Subroutine 
CSORT generates the following information: 


ND(J) number of dipole modes sharing segment J 

MO(J,K\ list of dipoles sharing segment J 

NCM size of the compressed open-circuit impedance matrix 
N number of dipole modes of the complete systems 


IDSEG(J) the plate (number) that acts as the ground plane 
for segment J. 


The following quantities must be specified in the main program. 


ICC dimension the compressed matric C(I,J) 

ICJ dimension related to the number of dipole modes N 
INP dimension related to the number of points NP 

INS dimension related to the number of segments 

NOP number of plates of the structure 

NOC number of corners of the structure. 


The numerical values assigned to the above quantities must 
agree with the dimensions actually reserved for the corresponding 
quantities in the COMPLEX and DIMENSION statements. For example 
when the structure has only one plate NOP=1, correspondingly we 
must redimension all related vectors, é.g., 

DIMENSION NPLC(1,4), VNP(1,4) 


11 


SSS 


instead of 
DIMENSION NPLC(10,4), VNP(10,4) for the 10-plate case. 


X(I), Y(I), Z(1) denotes the dimensionless quantities kx, ky, 
kz for point I where k = 2n/A. XC(I), YC(I), ZC(1) denotes 
the coordinates x, y, z in meters. 


The main program calls subroutine CMATX (Appendix 10) to gen- 
erate the compressed open circuit impedance matrix C(I,J). CMATX 
performs the same function as subroutine IDANT of Reference 3. The 
main difference lies with the addition of coupling terms due to 
reflected and diffracted fields in the matrix. 


The program next calls subroutine ANTI (reference 3) to obtain 
the resultant current distribution on the wires and the radiation 
efficiency EFF. If the wire antennas have only one generator and 
VG = (1.,0.) then Yaa “1 denotes the antenna input admittance 
and impedance respectively (see example 1, section IV). If there 
are more than one generator; input impedance at input port J can 
be obtained by dividing VG(J) by CJ(J) (see example 2, section IV). 
(Note that if the generator excites a point on the ground plane 
- CJ(J) should be used, as the current for the ground plane dipole 
mode is defined in reverse direction from those of the real modes 
(see also reference 3)). 


Finally the antenna pattern is obtained by calling subroutine 
CIFFLD (appendix 11). TH and PH denote the spherical coordinates 
6, in degrees of the distant observer. CIFFLD is called once for 
each observed direction, TH, PH. CIFFLD is essentially the same 
as subroutine JFFLD of reference 3 except for the call CZFF. Pat- 
terns in xy, yx, xy planes are obtained by setting pattern indi- 
cators IXY, IYZ, IXY respectively to a positive integer value, if 
all indicators assume zero or negative integer values no pattern 
is computed. The resultant patterns are stored in library file 
PLDATB ready for plotting. 
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| Figure 4. Main program. 
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Figure 4 (Cont'd). 
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Figure 4 (Cont'd). 
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IV. EXAMPLES 


A. A Quarter Wave Monopole on Plate 1 


The monopole is located at (0.,-A,0.) on plate 1 of the struc- 
ture in Figure 1. It is fed with 1 volt delta gap generator at 
the point of contact with the ground plane. The output file DAd1 
aS generated by the program is listed in Figure 5. Input impedance 
244 is (40.12+j21.8) ohms. If the same monopole is located at the 
center of a square plate 200\ each side, computation with the present 
program yields atl = 39.66+j21.54 ohms, which is the same value 
as that of a quarter wave monopole mounted on an infinite ground 
plane obtained with programs of reference 3. Far field patterns 
in x-z, and y-z planes are shown in Figures 6 and 7, without in- 
clusion of doubly diffracted fields. These patterns agree well 
with those independently generated by R. Marhefka's GTD computer 


program’. 
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B. Two Quarter Wave Monopotes on Plate 1 and Plate ? 


Consider the second example where the first monopole is lo- 
cated as in example 1. The second monopole is at center of face 
? or at (-2d,0.,-0.5X) in x, y, z coordinate system as shown in 
Figure 1. The 2nd monopole is perpendicular to its ground plane 
and fed with 1.0 volt delta gap generator at ground plane end. Out- 
put file DApl resulted from this program is listed in Figure 8, 
where input impedances for monopole 1 and 2 are yy and Zoo respec- 
tively. 


More importantly, the knowledge of the coupling or interference 
effects between wire antennas mounted on the satellite is essential 
in the design stage to evaluate alternative arrangements. The present 
program is well equipped for this purpose. For example, let us 
consider the case of two quarter wave monopoles mounted on plates 
1 and ? of the satellite (Figure 9). Here we wish to study the 
amount of interference or coupling between the two monopoles. Mono- 
pole on plate 1 is excited with a1 volt delta-gap generator, and 


the program (DOANI}, 3468N) computes the complex voltage Vio which 


appears across the 50 termination of receiving monopole on plate 
2. & prot of Vio versus distance d, is shown in Figure 10, where 
7) denotes the separation of the receiving monopole from the edge. 


The dB plot shows a change in the slope of IVy oI close to d, 3 
1.08, where the receiving monopole starts seeing the direct incident 
field from the monopole on plate 1. 
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APPENDIX J 
SUBROUTINE CROSSP 
CALL statement 
CALL CROSSP (V1, V2, VN3, V3) 


This subroutine computes the vector of products of V1 and V2 


— — a 
V3 = V1 x V2 
—— 
az V3 : 
VN3 = —— unit vector of cross product. 


al 


The listing of CROSSP is shown in Figure A.1. 
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Figure Al. Subroutine CROSSP. 


APPENDIX 2 
SUBROUTINE HELIX 


This subroutine is listed in Figure A.2. It describes the 


geometry of helical antennas. 


CALL statement 


CALL HELIX (XC,YC,ZC,1A,1B,CIRCM,ALFA, TURN,NANT,AXIS,ASTAR,NRS,NRP) 


INPUT. 

CIRCM: Circumference of helix 
ALFA: Pitch angle 

TURN: Number of turns 

NANT: Number of helical antennas 


AIXS(NANT,3): Position vector of the feed point of the helical 
antenna NANT 

ASTAR(NANT,3): Position vector of the standing point of the helical 
antenna NANT (or end point of the connecting segment) 


OUTPUT 

XC,YC,ZC: Coordinates of end points of all segments of the 
helical antennas 

IA, IB: List of index of end points 


NRS: Number of real segments | 
NRP: Number of real points. 
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Figure A2. Subroutine HELIX. 
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APPENDIX 3 
SUBROUTINE COMED 


This subroutine is listed in Figure A.3. It is used to find 
the plate (number) that shares corners Cy. Cy with plate I. 
CALL statement 


CALL COMED (1,C1,C2,CR,NPLC,NCED,NALL NOP ,NOC) 


Input: 

I: plate (number) under consideration 

C162: position vectors of corners C1, C2. Dimension C1(3),C2f2) 

CR: List of position vectors of corners. Dimension CR(NOC,3), 
where the first parameter NOC is the corner index. 

NPLC: List of corner indices. Each corner index carries two 


parameter NPLC (NOP,J) where NOP is the plate index and 
J is the integer 1, 2, 3 or 4 assiqned to each corner 
of the 4-corner plate. 

NALL>O: check out all plates (NALL plates) 

NALL<O: check only the plates with index number smaller than I 

NOP: total number of plates. 


Output: 
NCED>0 plate having index NCED shares corners C120, with plate I 


NCED<O: no plate shares corners C1 2Co with plate I. 
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Figure A3. Subroutine COMED. 
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APPENDIX 4 
SUBROUTINE IMAGE 


This subroutine is listed in Figure A4. Subroutine IMAGE finds 
the position vector RI(3) of the image of a source vector RS(3) 
projected on the plate (number) I. 


CALL statement 


CALL IMAGE (RS,RI,1,NIMA,VNP,NPLC,CR,NOP,NOC) 


Input: 

1 plate under consideration 

RS: position vector of source point 
NP: list of unit vectors normal to plate I. 
NPLC: list of corner indices. 

CR: list of position vectors of corners 
NOP: total number of plates. 

NOC: total number of corners. 

Output: 

RI: position vector of image point 
NIMA: image indicator 


NIMA>0 image exists and distinct from source 
NIMA<O source is on the plate, RS=RI. 
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Figure A4. Subroutine IMAGE. 
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APPENDIX 5 
SUBROUTINE REGION 


This subroutine is listed in Figure A.5. Given a source point RRS 
and an ohservation point RRO, subroutine REGION will check if any plate 
is blocking the ray path. A plate does not block the ray path if either 

the source point or the observation point is on the plate. 


CALL statement: 


CALL RFGION (RRS,RRO,RRT,CR,VNP,NNI,NN?, TELL,NPLC,NFAR,NALL MM, 


Ea an i NA ni mh ie a AE aR ci a i 


NOP ,NOC) 
Input: Position vectors of source point and observation point 
respectively. 
CR: List of position vectors of corners 
VNP: List of unit vectors normal to the plates 


NN1,NN2: Blocking check of plates NN1, NN2 not required. If want 
to check blocking by all plates, set NNI, NN2 to negative 


integers. 

NPLC: List of corner indices 

NFAR: Far field indicator. For observation point in the far 
field NFAR>O 

NALL: Number of plates that may block the ray path (usually 
NALL=NOP ) 

NOP: Total number of plates on structure 

NOC: Total number of corners on structure 

Output: 

RRT(3): Position vector of point of intersection 

MM: The plate (number) that first blocks the ray. 

TELLS Blocking indicator 


TELL=.TRUE.: A plate is blocking 
TELL=.FALSE. No plate is blocking the ray. 
Subroutine called is CROSSP. 
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Figure A5. Subroutine REGION. 
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APPENDIX 6 
SUBROUTINE DIFPT 


This subroutine is listed in Figure A.6. Given a source point 
RS and a near-zone observation point RO, subroutine DIFPT finds 
the position of the diffraction point on the edge (wedae) following 
the procedure descrihed in reference 1. 


CALL statement: 


CALL OTFPT (RS,RO,C1,C2,VNN,SP,S,ANI,RD,VI,IDX,VPHNS,VPHND,VBNS, 
VBND,PHS.PHD,BN,E) 


Input: 

RS(3),R0(3): Position vectors of source point and near-field 
observation point respectively. 

C1(3),C2(3): Position vectors of corners C1, C2 defining the edge 


VNN(3): Unit vector normal to the plate, ni. 

Output: 

RD(3): Position vector of diffraction point Q 

SP: Distance between source and diffraction point in meters. 

Se Distance between observation point and diffraction 
point in meters. 

ANI: Angle of incidence w.r.t. the edge in radians 

VI(3); Unit vector alonq incident ray from source to dif- 
fraction point 

IDX: Edae diffraction indicator IDX<0 no diffraction by 


the edge, vice versa 


VPHNS (3), VPHND(3) 8,8 unit vectors (see accompanying Figure, Figure 11) 
VBNS(3),VBND(3): 4', 4 unit vectors 


PHS,PHD: b',d angles in degrees. 
a i 9 Fe Unit vector along edge: from corner C1 to C2, E. 
BN(3): BN=VNNxE, unit vector perpendicular to the edge, B. 


Subroutine called is CROSSP. 


| Figure 11. Unit vectors for ray-fixed coordinate system. 
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Figure A6. Subroutine DIFPT. 
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APPENDIX 7 
SUBROUTINE DIFPTF 


This subroutine is listed in Figure A.7. Given a source point 
RS and a far-zone unit vector RO, subroutine DIFPTF finds the posi- 
tion vector of the diffraction point RD, if any, on the edge defined 
hy corner Cl and C?. Relevant expressions are described in Ref- 
erence 1. 


CALL Statement: 


CALL DIFPTF (RS,C1,C2,VNN,RO,SP,ANI,RD, IDX, VPHNS,VPHNS,VBNS, VBND, 


PHS,PHD,BN,E) 

Input: 

RS(3),RO(3): Position vectors of source point and far-zone 
unit vector 

€1,(3') €2 (3): Position vectors of 2 corners defining the edge 

VNN(3): Unit vector normal to the plate 

Output: 

RD(3): Position vector of the diffraction point on 
the edge. 

SP: Distance between source and diffraction point 
in meters. 

ANI: Angle of incidence w.r.t. the edge in radians 

IDX: Edge diffraction indicator; IDX<0, no diffraction 


by the edge, vice versa. 
VPHNS(3),VPHND(3);  $',6 


Aa 


VBNS(3),VBND(3); | B',g 


PHS ,PHD: 6'56 

BN(3): BN=VNNxE; also a unit vector perpendicular to 
the edge. 

E(3): Unit vector directing along the edge from corner 
Cy to Co. 
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APPENDIX 8 
SUBROUTINE WEDIFF 


Subroutine WEDIFF is listed in Figure A-8. Given a source 
point RS and an observation point RO in the far-field (NFAR>0O) or 
near-field zone (NFAR<0), subroutine WEDIFF generates relevant pa- 
rameters for the computation of the diffracted field in the presence 
of the edge defined by corners C1 and C2. 


CALL statement: 


CALL WEDIFF (RS,C1,C2,RO,RD,IPL,NFAR,AM,WAVM,NCM,VNP,CR,NPLC,NOP, 
NOC ,VBNS,VPHNS,CONS1,VB2,VB2, IDIF,VBND,VPHND) 


Input: 
RS(3),RO(3): Position vectors of source point and observation 
point 
IPL: Plate under consideration 
NFAR: Far field indicator 
AM: Wire radius in meter 
WAVM: Wavelength in meter 
NCM: Size of the compressed matrix 
VNP: List of unit vectors normal to the plates 
CR: List of position vectors of corners 
NPLC: List of corner index 
NOP: Number of plates 
, NOC: Number of corners 
| C1,C2: Corners Cl and C2 defining the edge 
Output: 


VBNS(3),VPHNS(3): Unit vectors 8',6' (Figure 11) 
VBND(3),VPHND(2): Unit vectors B',6 (Figure 11) 


IDIF: Diffraction indicator 


CONS]: Diffraction coefficient ' 
VB? , VB? VallsB sn), Vp(L,B an); B°=6#0! 


such that 


EF) = -(VB1-vB2)*E,"*CONS1 
1 = ~(VB1+VB2)*E, *CONSI 


mm 
Q 


Phase of ef and p is referred to origin of coordinate system when 
NFAR>0, to source point when NFAR<O. Subroutines called DIFPT, 
REGION, DIFPTF,COMED,CROSSP, 


CeeOFUCE UIFFR ACTION 
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Figure A8. Subroutine WEDIFF. 
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APPENDIX 9 
SUBROUTINE CSORT 


This subroutine is listed in Figure A.9. It is used to set 
up the image points and image segments, generate ID indices for 
these points and seqments, derive the size of the compressed im- 
pedance matrix, set up the modes and provide a list of modes MD(J,K) 
and a list of dipoles ND(J) sharing segment J. Basically CSORT 
is a modified version of subroutine ISORT in reference 2 to deal 
with wire antennas on different plates of the satellite structure. 
Input data for the satellite structure are entered through the last 
seven parameters in the CALL statement, namely VNP,NPLC,CR,NOP,NOC, 
IDSEG,IDPT. Output parameters are NCM,IDSEG,MD,N,0C,I1,12,13,IA(K), 
IB(K),XC(I),YC(1I),ZC(1). 


CALL statement: 


CALL CSORT (IA,IB,ICC,ICJ, INS, IWRITE,12,12,13,JA,JB,MAX,MIN,MD,N, 
NCM,ND,NP,NS,NRP,NRS,NPGP ,NSGP,DC,XC,YC,ZC,VNP,NPLC,RC,NOP,NOC,IDSEG, 
IDPT) 


The above parameters are defined as follows: 


Variables 


Input: 

IA( J), IB(J): End points index of segment J extending from 
A to B 

ICC: Dimension relating to matrix C(I,J) 

ICU: Dimensions relating to the number of dipole 
modes N 

INP: Dimension relating to the number of points 
NP 

INS: Dimension relating to the number of segments 
NS 
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XG. ¥E-76- 
NRP: 
NRS: 
NPGP: 


NSGP: 


NP: 


NS: 


VNP ,NPLC,CR,NOP NOP: 


Output: 
NCM: 
ND): 
MD(J,K): 
N: 
T,,1,,1,: 
XG, YC ,263 
IA, IA: 
JA, JB: 
DC: 
IDSEG(K); 
IDPT(1): 


Coordinates of real points in meters 
Number of real points 

Number of real segments 

number of points touching the qround plane 
(plates) 

Number of segments with end points on the 
plates 

Number of points for the complete antenna 
system 

Number of segments for the complete antenna 
system 

Parameters relating to satellite structure 
as defined previously 


Size of the compressed matrix 

List of dipoles sharing segment J 

List of dipole modes sharing seaqment J 

Number of dipole modes of the complete system 
Terminal point I,, end points , I, of dipole 
modes 

Coordinates of imaginary and real points 

End points of real and imaginary segments 
Segment number for each dipole mode 

Segment length in meters 

Plate that accomodates segment K 

Plate that relates to endpoint I 


Subroutine called IMAGE. 


IA, IA: 


DC: 
IDSEG(K); 
IDPT(I): 


Subroutine called IMAGE. 


Size of the compressed matrix 

List of dipoles sharing segment J 

List of dipole mode sharing segment J 

Number of dipole modes of the complete system 
Terminal point Io, end points I,, I, of dipole 
modes 

Coordinates of imaginary and real points 

End points of real and imaginary segments 
Segment number for each dipole mode 

Segment Tength in meters 

Plate that accomodates segment K 

Plate that relates to endpoint I 
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Figure A9. Subroutine CSORT. 
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APPENDIX 10 
SUBROUTINE CMATX 


Subroutine CMATX is listed in Figure A.10. It computes C matrix 

for a general wire configuration similar to subroutine IDANT and 

7GS of reference 3, but CMATX also accounts for the components of 

the C matrix due to reflected and diffracted fields from plates 

and wedges as expressed in the P(I,J) terms. In the form listed, 
CMATX does not include effects of wire conductivity. However, if 
wire surface impedance is to be included, CMATX can be extended 

hy adding the section from statements 200 and 260 of subroutine 
IDANT? before statement 26? of CMATX. 


CALL statement: 


CALL CMATX (ICC,INT,JPP,MD,N,NCM,ND,NLF,NP,NPGP,NRS,NS,AK,CMM,D, 
FMC ,CDK,SDK,NFAR,WAVM,AM,XC,YC,ZC,CR,VNP,NOP,NOC,NSGP, 
NPLC,X,Y,Z, IDSEG) 


Input and output parameters are defined below. 


Input: 
ICC Dimension related to the open-circuit impedance matrix C 
INT: Integer denotes the number of integration intervals using 
Simpson's rule. For general purpose set INT=4 
JPP: Number of modes that have images 
MD: List of dipole modes sharing segment J 
N: Total number of dipole modes 
NCM: Dimension of the compressed C matrix 
ND: List of dipoles sharing segment J | 
NLD: Number of lump loads | 
NP: Number of points 
NPGP: Number of points on the ground plane 
NRS: Number of real segments 
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NS: Number of segments 

AK: (2m/d,)a where a iS wire radius 

CMM: Wire conductivity in megamhos/m 

D: List of dimensionless segment lengths kd 

FMC: Frequency in MHz 

NFAR: Far field indicator 

WAYM: Wavelenath in meters 

AM: Wire radius in meters 

XC,YC,7C: Coordinates of the end points 

CR: List of position vectors of corners in meters 

VNP: List of unit vector normal to the plates 

NOP: Number of plates 
NOC: Number of corners | 
NSGP: Number of segments on the ground plane 

NPLC: List of corner indices 


XC,YC,ZC: Dimensionless coordinates of end points 


IDSEG: List of plates (numbers) that act as a ground plane for 
segment J. 

Output: 

CDK: List of cos(yd) for all segments 

SDK: List of sin (yd) for all seqments 

C: Resultant C matrix 


Subroutine called: 7GS,7GMM, REGION, IMAGE,COMED,WEDIFF 


sh 


SUBROUTINE CMATX 


DO 200 


SELECT SOURCE SEGMENT K 
AND RECEIVING SEGMENT L 
AND APPROPIATE MODES I,J 


K AND L ON THE 
SAME PLATE ? 


YES NO 


BOTH K & L REAL 
SEGMENTS ? 
NO 


CHECK DIRECT RAY 


COMPUTE COUPLING 
DUE TO DIRECT RAY 


120 
K &L SHARE 
CEE 2am ONE POINT ? 
CALL CALL 
ZGMM ZGS 


K&L REAL 
SEGMENTS ? 


CHECK REFLECTION 


COMPUTE COUPLING 
DUE TO REFLECTION 


DO 270 
CALL GNF 
270 CONTINUE 


CHECK DIFFRACTION 
RAY 


COMPUTE COUPLING 
DUE TO DIFFRACTION 


DO 100 


17 COMPUTE CST 
168 COMPUTE C(I,J) 


200 CONTINUE 


RETURN 
END 
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25 ¢ COP CRIACHEL« (COKE OO.) 


Figure Al0. Subroutine CMATX. 
57 


be ef CONTINUE 


ams TECOMITE LT .AK)GO Te 71 
| era TEIOMAK COTA) GM TO e1 
34 TrAKe T2020 TH Al 
= Ot) tu TO 2? 
ol el VARITE Chee AR et MAK TM 
6, Woh 
4 RETURN 
64 2? 0 200 KF1,5S 
6s NDK=hiDCK) 
bb KA=IA(K) 
67 “BSIB(K) 
bt mK=0 CK) 
6% SGNSSEC PL X (Ue gSDK(K)) 
Ti CONSEC™PLACCUK CK) eC.) 
71 bu 200 LolefS 
Te MUL =NDC(L) 
e343 tASIAG) 
7h PurIKttd 
tS meeb(l) 
Te SONTSC"PLX (0. eSDKILD) 
77 N{Ls=b 
76 NG 200 Tizielnk 
79 T="O(KeTl) 
BN LECIestelice)Go TO 206 
o4 Fi=l. 
62 TF (KB ef ele(L))G0 TO 3¢€ 
84 TR(KReot Ge TL (1) IFLl=ale 
84 Ts=1 
B85 Gu To 6a 
66 36 TE (Ke UeJ StF) Flsel, 
87 1s=2 
Sb 41 PO 200 JutVetpl 
8% J=a"D(LeJdu) 
Fu JECT ete) FO EhE 
9) PueLs 
é TE (LRei GeTely))G0 TO HE 
93 TE (LR ef Gell (Jd) FUEo1. 
a4 Joel 
3S co TO “6 
Su 4Eé TE( LAS! Qe S(J) )FURHl, 
97 Js=2 
ee TEONILe NE e0)GO TC 16% 
39 NJL=L 
200 C¥EePEETIVLiNG SIMILAR TO 268 
101 CAHOX(E REX CL ADIZOL 
10¢ CBHSCY (LED HKYOLADIZOL 
103 CORIZUL PD HEZOLADIZOL 
194 TNSECRECINT SZ?) 
165 TE CINSeLT eZ) INSZ2 
10¢ PP=LiSt1 
107 NMELTSOL/ZINS 
104 T=0. 
109 PLIT=A(U)eqeV) 
11u bLP=(Oeeel) 


Figure A10 (Cont'd) 
58 


ON 


Rel=(UseU) 
Per=(UeeeU) 
bHLITE (User) 
WLP=H(Ueee GD) 
HELE (UeeeW) 
P2P=(Usee") 
OlleLI=(U.s.0) 
(Ve2d=(Ueeel) 
C(2eL)=(b,,.0) 
COPe2d=(Ueee D) 
SGlizeLle 


CeeL TET UP TYE SULRCFS CORRESPONDING Te ENDPOINTS TOUCHING THE 
CetGE OHNE PLONE FCF BPEFLECTION ati OLFFRACTION TESTS, 
CHEN Le FAC ACCOUNTS FOR GRAZING INCTOE ACE 


20k 


21ly 


215 


2?) 


Dn 


TE CKAeoT eM eANDe KBeGTSPEP) GO TO 210 
TE (KBeLE oF GPIGU TO 208 
TUSEGK=ALIUSEGCK) 
RSACLISXC(KAD +1 eb “N6*VNE(LUSEOK 61) 
PSA(2V=YC (KE) +Lebk -UG*eVNP(ITLUSFGH,?) 
USA(SIAZCCKS Det eof CEeVNP( LUSFGK,3) 
RSP (1I=KCCK! ) 

RSP (2I=YC (KE) 

LSC (3I=Z2C(KY) 

AFACZU,S 

PFACH1,. 

66 TO 41% 

TUSEGHK=INSECIK) 

PSACLI=ENXC UKE) 

PSA(2)=YC(KA) 

ROA CAI=ZE (KAD 

PSE CLI=XC COKE DHT ek MCEeVNP CL OUSEGK 1) 
PSEEPVAYC CRE FL ek me U6eVNP( TOSFGK,2) 
RSP(SIEZC (KE) $1 eb er O*VNP(IUSFGK. 3) 
AFACH1. 

PFACH0.S 

GO TO “15 

PSACLIA=KCCHS) 

PSAC2I=YCURA) 

PSA(3)=ZC (KA) 

PSF(1L)I=XC (KP) 

PSP (AI=YCOK®) 

FSP (S)S7CCkK%) 

ferCHle 

PFACH1,. 

CASH(X (KE) eX (KANDI SOK 
CESZCY (KB) Re Y(IKA)I/(IK 
COF=H(ZIKE eZ (KADISIK 

COCSCAKCASSC RACBS+CG*CES 
TONHLUSE COR DRT OSEG(L) 

TREIDD MES IEO TO e4b 

TE (Wer Meh dCU TO 124 

TWHECL AKA) * (LBOKB) *(LAMevi) e( LBA) 
TECIN eFE2G)G0 TO FO 


CEeANTE NA SEGMENTS ON SAME PLATE 


CALL Z5SOX (KAD ATOKA) Oe Z0KA) e¥OKB) eV(KR ID Oe 7(KB) e¥(LAIOVILAD © 
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| APPENDIX 11 
SUBROUTINE CIFFLD 


This subroutine is listed in Figure A.11 and is essentially 
the same as subroutine IFFELD of reference 2, except for the CALL 
CZFF statement, where CZFF calculates the far zone field of a mono- 
pole above the satellite structures (see Appendix A.11). 


Let (r,9,) denote the spherical coordinates of the distant 
observer, and let E (I), I (I) denote the electric field intensities 
of dipole mode I with unit current. Then 


EPP(I) =re 4 
ETE(I) = ee” ECT) 
Summing the fields due to all modes I yield 


N 
EPH = } CU(1)EPP(I) 


N 
ETH } CO(IT)ETT(I) 


where CJ(I) denote the terminal current of mode I 
CALL statement 


CALL CIFFLD (INS,MD,N,ND,NRS,CDK,CJ,DC,EPH,ETH,G,CPP,GTT,PH, SDK, 
TH,XC,YC,ZC,WAVM,AM,NPGP ,NCM,VNP,NPLC,CR,NOP,NOC, IDSEG) 


All symbols have been defined in previous subroutines. 


A5 
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APPENDIX 12 
SUBROUTINE CZFF 


CZFF is listed in Figure A.12. 


Submatrix CZFF calculates the far zone field of a monopole 
above a satellite body modeled by flat plates. CZFF includes con- 
tributions from incident, reflected, singly diffracted and doubly 
diffracted fields. All phase terms of these fields are referred 
to the origin of the coordinate system before summation. 


Between statement of 26 and 30 the direct incident field region 

is defined and the incident field computed by calling subroutine 

GFF (reference 3). The fields due to reflected rays are computed 

in loop Do 5. Do 6 defines diffraction reqion and computer singly 
diffracted fields. Doubly diffracted fields are evaluated between 
statements 80 and 40. Finally the total field arrives at the ob- 
servation point due to each end of the monopole is the sum of the 
above field components. Summation is performed between statement 

40 and RETURN. 
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ET1,EP1 are 6 and » components of the far field due to endpoint 
1, E72, EP2 and those due to endpoint 2. Other symbols have been 
different previously. 
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